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Abstract 

The expansion of space, and other geometric properties of cosmological models, can 
be studied using geometrically defined notions of relative velocity. In this paper, we con- 
sider test particles undergoing radial motion relative to comoving (geodesic) observers in 
Robertson- Walker cosmologies, whose scale factors are increasing functions of cosmolog- 
ical time. Analytical and numerical comparisons of the Fermi, kinematic, astrometric, 
and the spectroscopic relative velocities of test particles are given under general circum- 
stances. Examples include recessional comoving test particles in the de Sitter universe, 
the radiation-dominated universe, and the matter-dominated universe. Three distinct 
coordinate charts, each with different notions of simultaneity, are employed in the cal- 
culations. It is shown that the astrometric relative velocity of a radially receding test 
particle cannot be superluminal in any expanding Robertson- Walker spacetime. How- 
ever, necessary and sufficient conditions are given for the existence of superluminal Fermi 
speeds, and it is shown how the four concepts of relative velocity determine geometric 
properties of the spacetime. 

1 Introduction 

General relativity provides no a priori definition of relative velocities of non local objects 
in curved spacetime. Different coordinate systems and notions of relative velocity yield 
different results for the motion of distant test particles relative to a particular observer. This 
ambiguity led to consideration of the need for a strict definition of "radial velocity" at the 
General Assembly of the International Astronomical Union (IAU), held in 2000 (see [TJ [2]). 

Thereafter, a series of papers [3j|4j|5j appeared addressing the general question of relative 
velocities and culminated in the introduction of four geometrically defined (but incquivalent) 
notions of relative velocity: Fermi, kinematic, astrometric, and the spectroscopic relative 
velocities. These four relative velocities each have physical justifications, and have been 
employed for studying properties of spacetimes [HI |Tj . Related work includes [5] for the study 
of observer-referred kinematics and dynamics, [10 for a study of Schwarzschild black holes, 
and [HJ [12] for Doppler tracking. 

The four definitions of relative velocities depend on two different notions of simultaneity: 
"spacelike simultaneity" (or "Fermi simultaneity" [13]) as defined by Fermi coordinates of 
the observer, and "lightlike simultaneity" as defined by optical (or observational) coordinates 
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of the observer [14] . The Fermi and kinematic relative velocities can be described in terms 
of the former, according to which events are simultaneous if they lie on the same space 
slice determined by Fermi coordinates. The kinematic relative velocity is found by first 
parallel transporting the 4- velocity v! of the test particle at the spacetime point q s , along a 
radial spacelike geodesic (lying on a Fermi space slice) to a 4-velocity denoted by T qaP u' in 
the tangent space of the central observer at spacetime point p, whose 4-velocity is u. The 
kinematic relative velocity i>km is then the unique vector orthogonal to u, in the tangent 
space of the observer, satisfying T qsP u' = j(u + i>kin) for some scalar 7 (which is uniquely 
determined). For a test particle undergoing radial motion, the Fermi relative velocity, VFermi 
is the rate of change of proper distance of the test particle away from the central observer 
along the Fermi space slice, with respect to proper time of the observer. 

The spectroscopic (or barycentric) and astrometric relative velocities can be found from 
spectroscopic and astronomical observations. Mathematically, both rely on the notion of 
"lightlike simultaneity", according to which two events are simultaneous if they both lie 
past-pointing horismos (which is tangent to the backward light cone) at the spacetime point 
p of the central observer. The spectroscopic relative velocity u sp ec is calculated analogously 
to Wkin, described in the preceding paragraph, except that the 4-velocity u' of the test particle 
is parallel transported to the tangent space of the observer along a null geodesic lying on the 
past-pointing horismos of the observer, instead of along the Fermi space slice. The astrometric 
relative velocity, u as t , of a test particle whose motion is purely radial is calculated analogously 
to i>Fermij as the rate of change of the affine distance, which corresponds to the observed proper 
distance (through light signals at the time of observation) with respect to the proper time 
of the observer, as may be done via parallax measurements. We describe this more precisely 
in the sequel, and complete definitions for arbitrary (not necessarily radial) motion may be 
found in [5]. 

In [7j exact Fermi coordinates were found for expanding Robertson- Walker spacetimes and 
were shown to be global in the non inflationary case. Fermi coordinates were then used to 
calculate the (finite) diameter of the Fermi space slice, as a function of the observer's proper 
time, and Fermi velocities of (receding) comoving test particles. In this paper, we extend the 
results of [7] by calculating all four relative velocities for test particles undergoing arbitrary 
radial motion in general expanding Robertson- Walker spacetimes. Examples include explicit 
expressions for the Fermi, kinematic, astrometric, and the spectroscopic relative velocities 
of comoving test particles in the de Sitter universe, the radiation-dominated universe, the 
matter-dominated universe, and more generally, cosmologies for which the scale factor, a(t) = 
t a with < a < 1 (see equation ([I]) below). 

We express the metric tensor in optical coordinates to calculate the two relative velocities 
associated with lightlike simultaneity. Fermi coordinates are utilized to find functional rela- 
tionships between all four relative velocities. Necessary and sufficient conditions are given 
for the existence of superluminal Fermi relative velocities of test particles undergoing radial 
motion, and the astrometric relative velocity of a radially receding test particle is shown to 
be necessarily subluminal (the other two relative velocities are subluminal by their defini- 
tions). In addition we show how the relative velocities determine the leading coefficients of 
the metric tensor in Fermi and optical coordinates (denoted as g TT and g TT below) and that 
pairs of them determine the scale factor, a(t), of the Robertson- Walker spacetime. 

This paper is organized as follows. In Section [2j we establish notation, define the four 
relative velocities for our circumstances, express the metric in Fermi and optical coordinate 
systems, derive general properties of the relative velocities and explain how they may be 
compared. In Section [3] we specialize to the case of test particles comoving with the Hubble 
flow and then apply the formulas in Section [4] to the Milne universe, the de Sitter universe, 
and to Robertson- Walker spacetimes for which the scale factor has the form a(t) = t a with 
< a < 1. This latter class of cosmologies includes the radiation-dominated and matter- 
dominated universes. In Section [5j we again consider general Robertson- Walker spacetimes 
and find functional relationships between the four relative velocities and the underlying ge- 
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ometry through the metric tensor. Section [6] gives concluding remarks with a discussion of 
relative velocities and expansion of space. 

2 Relative velocities 

The Robertson- Walker metric in curvature-normalized coordinates (or Robertson- Walker co- 
ordinates) is given by the line element 

d,s 2 = -dt 2 + a 2 (t) (d x 2 + sl ( X ) dn 2 ) , (1) 

where dfi = d9 + sin 2 Odip, a(t) is a positive and increasing scale factor, with t > 0, and 

(sin (x) if k = 1 
X if k = 

sinh (x) if fc = — 1. 

There is a coordinate singularity in at x = 0, but this will not affect the calculations 
that follow. Since our purpose is to study radial motion with respect to a central observer, 
it suffices to consider the 2-dimensional Robertson- Walker metric given by 

ds 2 = -dt 2 + a 2 (t)d X 2 , (2) 
for which there is no singularity at % = 0. 

2.1 Notation 

We denote a central observer by /3(r) := (r, 0), and a test particle by (3' (V) = (t (r') , x {?'))■ 
Both spacetime paths are parameterized by their proper times, and we will always assume 
that x > 0- O ur aim is to study the relative velocities of (3' with respect to, and observed 
by, (3. The 4-velocity of f3 is denoted by U := = (1, 0), and the 4-velocity of f3' is given by 
U' := ijj + X~§^ = (iy x): where the overdot represents differentiation with respect to r', the 
proper time of j3' . From g(U' , U') = —1, we obtain 

i = v/a 2 (i)x 2 + 1. (3) 

Vector fields will be represented by upper case letters, and vectors by lower case letters. 
Following this notation, the 4-velocity of (3 at a fixed event p = (t, 0) will be denoted by 
u = (1,0). The Fermi, kinematic, spectroscopic, and astrometric relative velocity vector 
fields (to be defined below) are denoted respectively as Vkin, ^Fcrmi, V spec and Vast- They are 
vector fields defined on the spacetime path, f3, i.e., the central observer. Since all of these 
relative velocities are spacelike and orthogonal to U, they are each proportional to the unit 
vector field S := ^F] 

Comparisons of the four relative velocities for a given test particle are possible. Direct 
comparisons may be made of the Fermi and kinematic relative velocities, because of the com- 
mon dependence of these two notions of relative velocity on spacelike simultaneity. Similarly, 
direct comparisons of the astrometric and spectroscopic relative velocities are also possible. 
However, a comparison of all four relative velocities made at a particular instant by the 
central observer (3 is possible only with data from two different spacetime events (q s and qe 
in Figure [lj of the test particle. Such a comparison, to which we refer as an instant com- 
parison, therefore lacks physical significance, unless the evolution of the test particle (3' can 
be deduced from its 4-velocity at one spacetime point, e.g. for comoving or, more generally, 
geodesic test particles. 

"This should not to be confused with the relative position vector field S used in [5]; in fact, S is the 
normalized version of S. 
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Figure 1: Scheme of the elements involved in the study of the relative velocities of a test 
particle (3' with respect to the central observer /3. The curves ip and ip* are spacelike geodesies 
orthogonal to the 4-velocity of /3, and A is a lightlike geodesic. 



It is also possible to compare all four relative velocities at a fixed spacetime event qi of the 
test particle through observations from two different times of the central observer, identified 
as r and t* in Figure [l] In the sequel, we refer to such a comparison of the relative velocities 
as a retarded comparison. 

In all that follows, we use the following notation for vectors at a given spacetime point 
p = (t, 0) and a spacetime point p* = (r*,0) in the past of p (see Figure [l]): 

^Fcrmi • — ^Fcrmi pi ^kin ' — ^kin p* j ^Fcrmi ' — ^Ferari p* 5 ^spcc • ^spcc p and ^ as t • — Vast p- 

So, in an instant comparison we compare Vkm, ^Fcrmi, w spcc , and w as t, while in a retarded 

and i; ast . 

2.2 Spacelike simultaneity and Fermi coordinates 

Let p = (r, 0) be an event of the central observer (3 with 4-velocity u. An event q is spacelike 
(or Fermi) simultaneous with p if ^(exp" 1 q,u) = fj^J The Fermi space slice Ai T (see [7], it 
is called Landau submanifold and denoted by L pu m [3]) is composed of all the events that 
are spacelike simultaneous with p. 
More explicitly, the vector field 



is geodesic, spacelike, unit, and X p is orthogonal to the 4-velocity u — (1,0) at p, i.e. X p 
is tangent to M. T . Let q s := (t s ,Xs) be the unique event of /?' D A4 T . Then there exists an 
integral curve of X from p to q s (the geodesic ip in Figure [T]), and so, using Q we can find 
a relationship between t, i s and 

a ( T ) . I(a(r)\ 2 _ 1 it s a(r) ^(rj-a 2 ^) 



From ([5]) we obtain that t s < r. Moreover, since we only consider positive coordinate times, 
we have to impose t s > and then it is necessary that < Xs < Xsmax, where 

r*<f>- 1 if. <6) 

TThe exponential map, exp p v, denotes the evaluation at affine parameter 1 of the geodesic starting at the 
point p, with initial derivative v. 
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On the other hand, with respect to time, it is necessary that t > r sm i n , where T sm i n is defined 
implicitly by 

a ( T ) 1 i- r~\ 

o a{r) ^a 2 (T)-a 2 (T) 

Note that x sm ax(T sm in(Xs)) = Xs and r smin (x sm a X (T)) = r. 

It will also be useful in what follows to consider Fermi coordinates. In [7J the metric 
([!]) was expressed in Fermi coordinates with respect to the central observer /3 (called Fermi 
observer in this context) for expanding Robertson- Walker spacetimes, and it was shown 
that Fermi coordinates are global when the spacetime is non inflationary. In two spacetime 
dimensions, the metric ([2]) expressed in Fermi coordinates (r, p) is given by 

ds 2 = g TT dT 2 + dp 2 , (8) 

where p is proper distance along a spacelike geodesic, orthogonal to the path /?, i.e., the 
Fermi distance from j3 to the corresponding event (see |5J). A formula for g TT in terms of 
(t, p) is given in [7j Theorem 2] . 

The vector field X given in Q may be expressed in Fermi coordinates as X — d/dp. By 
design of Fermi coordinates, r = t on the path /3(t) of the Fermi observer (where p = 0), 
but the two time coordinates differ away from that path. The coordinate transformations 
between Fermi coordinates (r, p) and curvature-normalized coordinates (t, \) & r e given as 
integral expressions in [JJ, but will not be needed here. 

Remark 2.1 Setting ds 2 = in ([8]), shows that the velocity of a distant photon with space- 
time coordinates (r,p), relative to the Fermi observer (3, is given by \dp/dr \ = \/—g T T{T, p). 
Thus, the metric ^ may be understood as a natural generalization of the Minkowski metric 
when the speed of a photon depends on its spacetime coordinates, and may expressed as 
ds 2 = — c(r, p) 2 dr 2 + dp 2 where c(r, p) is the Fermi speed of a photon at the spacetime point 
(r,p) relative to the Fermi observer located at (r, 0). 

2.3 Fermi and kinematic relative velocities 

In this section, we find general expressions for the Fermi and kinematic velocities of a radially 
moving test particle, relative to the central observer (3. These formulas will be applied to 
more specialized circumstances in later sections. 

Using the Fermi coordinates developed in the previous section, let u' s = rj^ | + p-j^ 

(f , p) be the 4- velocity of a radially moving test particle j3' at a point q s with Fermi coordinates 
(t, p), where the overdot represents differentiation with respect to proper time of (3', and 
p = (t, 0) is the event of from which we measure the velocities, with 4-velocity u = (1,0). 
Taking into account [51 Proposition 3], the Fermi relative velocity of u' s with respect to u is 
given by 

^Fermi = "T^p = V<5pj (9) 

dr t 

where S p = d/dp\ p . The requirement that g(u' s ,u' s ) — —1 forces ||wFormi|| < \J — 9tt(t, p), 
which, by Remark |2.1| has the physical interpretation that the right hand side of the in- 
equality is the speed, \dp/dr\, of a distant photon with spacetime coordinates (r,p), relative 
to the Fermi observer j3. It is an upper bound and limiting value for the Fermi speed of a 
massive particle. We therefore have: 

Proposition 2.1 In an expanding Roberts on- Walker spacetime, the Fermi relative velocity 
of a radially moving test particle at position (r, p) satisfies ||fFcrmi| < y-~9rr(j~, p) and can 
therefore exceed the central observer's local speed of light (c = 1) within a Fermi coordinate 
chart if and only if —g TT (T, p) > 1. 
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The following examples illustrate the proposition. 

Example 2.1 a) For Milne spacetime, — 3tt(t, p) = 1, the Fermi chart is global, and thus 
all Fermi relative speeds are subluminal. 

b) For the de Sitter universe, —g TT (r,p) — cos 2 (H p), with H p < n/2 (see [T51 UM [7]) 
where Hq is the Hubble constant. The Fermi chart is valid up to the cosmological 
horizon of this spacetime. Thus, all Fermi relative velocities are less than the local 
speed of light. 

c) For the radiation-dominated universe, i.e., for the case that a(t) — \fi in ([lj, the Fermi 
chart is global and 

1 2 

- 9rr (r, p) = - (1 + V<r - 1 sec" 1 ^a) , (10) 

where a > 1 depends on p and r. It may be shown that for any t > 0, the least upper 
bound of p is |r and that g Tr (r, p) — > | asymptotically as p — > |r 7 . Thus, 
Fermi relative speeds can exceed the speed of light in this spacetime, but are bounded 
above by |. 

An alternative expression for the Fermi relative velocity follows by expressing the 4- 
velocity of the test particle as 



(T,P) V-9rr(r,p) dr 



op 



(11) 



where ± indicates the sign of p, and £ := -J '—g TT (T, p) f is the energy per unit mass of the 
test particle as measured by an observer with fixed Fermi spatial coordinate p, i.e., by an 
observer with 4-velocity 

y/-g T r(T,p) dr ^ 

so that — g(uo,u' s ) = £. Note that a geodesic path for the test particle is determined by the 
value of £ at a single spacetime point. With this notation, we can write 

II^Fcrmill = V-9tt(t,p) ^ - . (13) 

Let Tq sP represent the parallel transport from q s to p by the unique geodesic joining q s 
and p, i.e. i(j in Figure [l] Using the relations g (T qsP u' s , d/dp\ p ) = g (u' s , d/dp\ qs ) and 
g (T qts pU' s ,T qa pu' s ) = —1, , we can find T qsP u' s , and then, by [3 Equation (4)], obtain the 
kinematic relative velocity of u' s with respect to u as 

v kin = —S p . (14) 

V-9tt(t, p) dr 

Thus, taking into account ^ and ( 13 1 we have 

hkinW = £ • (15) 

Returning to the curvature-normalized coordinates, and exploiting the symmetry of this 
spacetime through the killing field d/dx, the kinematic relative velocity may be expressed 
explicitly in terms of x s . From g (t w u' s , X p ) = g (u' s ,X qs ) and g(T qsP u' s ,T qs pu' s ) = -1 we can 
obtain T qsP u' s , and then find 

«kin = , = 'Op, (16) 



y(ts V /^S- 1 + a (^) 2 + 1 
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where t s = ^a 2 (t s ) X 2 s +l by (§. 

Also, from (|9| , the Fermi relative velocity of u' s with respect to u is given by 

dp f _i_ dp ,-, 

"Fermi — .<->p — g T ; . °p: V 1 1 ; 

where the function r(t s ,Xs) is defined implicitly by ((5j), and 

p(t s , Xs ) = dl°™(p, a,) = / 1 == df , (18) 

taking into account Q and [5J Proposition 2]. 
2.4 Lightlike simultaneity and optical coordinates 

Let p — (t, 0) be an event of the central observer (3. An event is lightlike simultaneous with 
p if it lies on the past-pointing horismos E~ (which is tangent to the backward light cone at 
the spacetime point p). 
The vector field 

_«(r)5_ o(r[_9 f . 

o(t) 9< a 2 (i) % v ; 

is geodesic, lightlike, and the integral curve A such that A(0) = p is a past-pointing null 
geodesic, affinely parameterized that satisfies the hypotheses of Proposition 6]. Let qt := 
[tii Xt) be the unique event of j3' n E~ . Then A is the unique geodesic from p to qg, and so, 
using (19) we can find a relationship between r, ^ and xt'- 

-1 Z" 7 " 1 

■df = / dx =► / — rdf = X l- (20) 



J J t a(r) 



From (20) it follows that ^ < r. Moreover, since we consider only positive coordinate times, 
tg > 0, and then it is necessary that < xt < X&nax( r )> where 

Xfin«x(r) := / - 1 - df (21) 

is the particle horizon for the observer /? at p. On the other hand, with respect to time, it is 
necessary that r > Tg m ; n (x^)) where T£ m - m (xe) is defined implicitly by 

1"fmin(X<!) 1 

^df = X , (22) 

Observe that Xtmax{rtmin(xt)) = Xt and r^ m i n (x^ max (r)) = r, so that each function is the 
inverse of the other. 

In the framework of lightlike simultaneity, it will be convenient to use optical (or obser- 
vational) coordinates with respect to the observer /3. Referring to Figure [I] we denote the 
optical coordinates of the point qg = (tg,xt) by (t,5), where the affine distance 8 from p to 
qi is defined as the norm of the projection of exp" 1 ^) onto the orthogonal complement u 1 - 
of u (see [5]). From (20), r(tg,xi) is determined implicitly, and differentiation gives 

dr cl(t) 8t , . . . 

dt e a(t e ) dxi 



It follows from ( 19 ) and [5J Proposition 6] that 

Xi) 

a{r{te,xi)) 



T(t t ,Xt) a ( f \ 

8 = 8{t t ,xt)= I , g ' » df ■ ( 24 ) 
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Differentiating ( 24 1 and using ( 23 ) , gives 



a(r) a(t e ) d(r) 
— o- 



06 



dig a{ti) o,(ti) 



85 



= a(r) — 5d(r), 



(25) 



where a{t) is the derivative of a(t). Now using (23 1 and (25 1, we may express the Robertson- 
Walker metric in optical coordinates (with respect to /3) in the form 

+ (26) 

a(T) 2 a z (T) J 



ds z = g TT dr z + 2drd5 = -2 1 



where te(r, 5) is given implicitly by (24) 



2.5 Astrometric and spectroscopic relative velocities 



Analogous to Section [273| we find here general expressions for the astrometric and spectro- 
scopic velocities of a radially moving test particle, relative to the central observer (3. The 
formulas will be used in later sections. 



Using the optical coordinates developed in the previous section, let u' e = r 



dr W 



' dSW 



(f , 5) be the 4-velocity of a radially moving test particle /3' at a point qe with optical coor- 
dinates (t, p) in the past-pointing horismos E~ , where the overdot represents differentiation 
with respect to proper time of j3' , and p = (r, 0) is the event of (3 from which we measure the 
velocities. From 5, Proposition 7], (26), and taking into account that g(u' e ,u' e ) = — 1, the 



astrometric relative velocity of u' e with respect to u is given by 



_dS _8 _ 1 
dr t 2 



1 



2 I S P- 



(27) 



There is no upper bound for ||w a st|| in the case of a radially approaching test particle (i.e. 



for the case dd/dr < 0) because f can be chosen to be arbitrarily close to zero in (27). 



However, a sharp upper bound for the case of a radially receding test particle is given by 
ll^astll < ~ 5tt/2, where the right side of this inequality is the relative speed, dS/dr, of a 
distant radially receding photon. Since it follows from (26) that — g TT < 2 when d(r) > 0, 
we have the following general result. 

Proposition 2.2 In any expanding Robertson- Walker spacetime, the astrometric relative 
velocity of a radially receding test particle is always less than the central observer's local 
speed of light (c = I). 



Using (26), it is easy to show that 6 is an affine parameter for the null geodesic segment, 
\(5) = (t, So — 5) for a given 5o > 0, where r is fixed and < 5 < So (see Figure[T]). Moreover, 
the vector field (19) is given by Y ~ d/88 in optical coordinates. Observe that by design of 
optical coordinates, f = along X(S). 

Let v, v' be the frequencies observed by u, u' e , respectively, of a photon emitted from the 
spacetime point q£. Then, the redshift of the test particle with radial motion relative to the 
central observer at p is determined by the frequency ratio 



v 

V 



(28) 



where P 
yields 



and thus, 



-d/dS is the 4-momentum tangent vector field of the emitted photon. Using (26) 



(29) 



(30) 



8 



We may also use (29) along with (26) in (27) to write an alternative expression for the 
astrometric relative velocity, 

«ast = 1 



6(t) 1 a 2 (t e ) 
a(r) 2 I a 2 (r) 



(31) 



In order to find expressions for the astrometric and spectroscopic relative velocities in 



terms of curvature-normalized coordinates, we make use of (|23| to obtain 

dr . 



v 
v 



dt 



(32) 



taking into account that if. — y 1 o?(ii)x\ + 1 by Combining this last expression with 
Q yields 



a 2 (T)(Vxi+a-Ht t )+ Xe ) -1 



Similarly, combining p2[) with (|27[) gives 



Wast, = 1 



&(t) a 2 (^) 



o(t) 2a 2 (r) 



1 + f^/o a (ti)x? + l + a(<£)x^ 



(33) 



(34) 



Remark 2.2 Expression (34) can be derived beginning with an expression analogous to (17) 
and using (23) and (25). 



Remark 2.3 For "retarded comparisons", the kinematic and Fermi relative velocities of u' e 
must be calculated relative to u* = Jj| » = (1,0), i.e., the 4-velocity of (3 at p* — (r*,0). 
From ([5]), r* = r* (te(T,Xi),Xt.) is defined implicitly from 



a(r*) 



dr = Xt, 



(35) 



where ti(T,xi) is given implicitly by (20). 



3 Comoving test particles 

In this section we apply the general results obtained in Section [2] to the case of test particles 
that are comoving with the Hubble flow. A comoving test particle parameterized by its proper 
time is given by /3'(r') = (r',x), where x > is constant, so that its 4-velocity U' = (1,0) 
and therefore u' s = (1,0) and u' £ — (1,0). Referring to Figure [lj we see that q s = (t s ,x) an d 
Qe = (te,x)- 

In ([5]) < s is implicitly defined as a function of (r, %), and similarly in (20) tf is implicitly 



defined as a function of (t, x)- In this section and the next, it will be convenient to regard 
not only t s and tg as functions of (r, x), but also the four relative velocities. However, it 
is important to recognize that in this context % is a parameter that labels a comoving test 
particle (with fixed coordinate x)i an d t is the time of observation by the central observer /3. 
The relative velocities are vectors in the tangent space of the point p = (r, 0) for test particles 
with coordinates (i s ,x) i n the case of the Fermi and kinematic relative velocities, and with 
coordinates (ti,x) m the case of the astrometric and spectroscopic relative velocities. Since 
all the velocities are proportional to S p and in the same direction, we will find expressions 
only for the moduli of the relative velocities. 

Remark 3.1 For the purpose of making retarded comparisons of comoving test particles, 



observe that t s (t*,x) — ti (r, x) follows from Remark 2.3 
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3.1 Kinematic and Fermi relative velocities of comoving test parti- 
cles 



Applying ( 16 ) with constant x, we obtain 

1 1 ^ kin 1 1 = 



a 2 (t s ) 



(36) 



Adapting (181 to the notation of this section, the Fermi distance from p = (t, 0) to 
9s = (t s (T, x), x) mav be expressed as a function of (t, x) as 



P = p( t , X) 



a(f) 



ts(T, X ) V a2 ( T ) ~ fl2 ( T ) 



df. 



Then, from 



\\VFc- 



dT k(r, x ) y/a*{r)-a*(f) 



(37) 



(38) 



3.2 Spectroscopic and astrometric relative velocities of comoving 
test particles 

Applying ( [33] ) with constant x gives 

a 2 ( T )-a 2 {t e ) 



IKpe 



Similarly, with x = 0, (34) becomes 

IKt|| = l-<5 



a 2 (r) + a 2 (t £ )- 

a(r) _ a 2 (U) 
a(r) a 2 (r) ' 



(39) 



(40) 



where 8 is the affine distance from p = (r, 0) to qg = (te(T,x),X)i anc i ( see (24)) it can be 
expressed as a function of (r, x)i 



6 = 8{t,x) 



a(f) 
a(r) 



df. 



(41) 



The astrometric speed has alternative formulations. The cosmological redshift of a co- 
moving test particle with respect to a central observer is well known and also follows from 



i/ air) 
v a(te) 



(42) 



Combining (39) and (42 1 in (40), we have the following result analogous to Hubble's Law for 
the Hubble speed. 

Proposition 3.1 For a test particle comoving with the Hubble flow and with redshift z, the 
astrometric relative speed relative to the central observer is given by 



I Vast I 



1-HS ' " Wspoc| 



1 + IKpecll 



1-HS 



1 



where || v spcc || is the spectroscopic relative speed and H = a(r)/a(r) is the Hubble parameter. 



10 



4 Comoving test particles in particular spacetimes 



In this section, we find explicit expressions for the relative velocities of comoving test particles 
for particular Robertson- Walker spacetimes. We consider the Milne universe, the de Sitter 
universe, and Robertson- Walker cosmologies with scale factors following a power law, i.e., 
a(t) = t a , where < a < 1, including, in particular, the radiation-dominated universe and 
the matter-dominated universe. As in the previous section, we continue to express t s and 
tg and the moduli of the four relative velocities as functions of (t, x). In fact for all of the 
examples discussed in this section, this dependence is exclusively through a parameter v given 

by 

v = v(t,x) := a(r)x, (43) 

where the overdot represents differentiation with respect to r. However, this is not true in 
general, e.g. if a(t) = t 2 +t. 



Remark 4.1 The expression (43) for v is the Hubble speed of a comoving test particle 



with curvature- normalized coordinates (r, x)- However, it is important to recognize that the 
relative velocities that we calculate in this section, as functions of v, are those of test particles 
located at different spacetime points, (t s ,x) and (tg,x)- 



Recall that for "retarded comparisons," defined in Section [2T| it is necessary to calculate 
the kinematic and Fermi relative velocities at p* — (r*, 0), which will be expressed in terms 
of the parameters (r*,x) instead of (r, x)- For this purpose, we define v* := v(t*,x)- 



4.1 The Milne universe 

The Milne universe may be identified as the forward light cone in Minkowski spacetime, 
foliated by negatively curved hyperboloids orthogonal to the time axis. As the four relative 
velocities in Minkowski spacetime were previously found in [5], we include this example only 
for purposes of illustration of techniques presented in this paper. 



For this spacetime, a(t) = t and k = —1, but our results are valid for any k. From (43) 

v = v(t,x)=X, ( 44 ) 

which does not depend on r. 



4.1.1 Spacelike simultaneity in the Milne universe 

By (|6]), XsmaxM = Jq ■= ^ T 2_ f 2 — +oo, and hence, by (44), v has no upper bound in the 
framework of spacelike simultaneity. From ^ 

" T 1 



T \It 2 — T- 



dr = x 



t s (r,x) 



coshx 



and so by (37) 



ts(r,x) Vt 2 - T 2 



dr = r tanh \ < r. 



From (38) 



dp ., p 

— = tanhx = tanhu = v k in = -■ 

OT T 



Applying (36) gives 




(45) 

(46) 
(47) 

(48) 



Thus, ||wkin|| = IkFermiH < 1- 
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4.1.2 Lightlike simultaneity in the Milne universe 



By (21), X£max( T ) — Jo h dr = +00, and hence, by (44), we have that v has no upper bound 
in the framework of lightlikc simultaneity. By ( 20 ) 



- dr = x 

. T 



So, it follows from (41) that 



S(t,x) 



— dr = re x sinh x- 

te(r,x) T 



Then, from (40) 



— — = e x sinh v = e v sinh v = —. 

OT T 



Similarly, (39) gives 

|| ^spec I 



T - tf 
T 2 +t 2 



= tanh x = tanh v = 



S/i 



I -5/i 



(49) 

(50) 
(51) 
(52) 



and it follows fro m |50l th at 8/t € [0, 1/2). 

From Remark 3.1. (45 1 and (|49|, we have r* = re~ x coshx, but for the Milne universe, a 



retarded comparison is equivalent to an instant comparison, because v* = v, and the relative 
velocities do not depend on r. 

4.2 The de Sitter universe 

For this spacetime, a(t) — e Hot , where Hq > 0, and k = 0, but our results are valid for any 
k. We suppose that t > 0, but we could also consider non-positive values for t. From (|43|) 

v = v(T,x)=H e H ° r X . 
4.2.1 Spacelike simultaneity for the de Sitter universe 

By ®, X < Xsma X (T) = j Q T -j^^=L= df = i-v/l - e - 2H « T for any point (i, x) on 

Fermi space slice at proper time r of the central observer. Thus, by (53), < v(t, ■) < 
^J e 2H T _ I and v ^ x ) > Hox 



(53) 



the 



By (§ 



,i?o(r-f) 



dr = x 



i 8 (r,x)=r-— ln(l + u 2 
^0 



and thus from (37) 



.Hot 



df 



Ho 



vT 



Then, from (38) 



Applying (36) and (54) gives 



I kid, 



dp v 
dr 1 + v 2 



o Hot 



sin (H p) cos [H p) ■ 



sin (if p) , 



VT+v 2 



(54) 

(55) 
(56) 

(57) 



where H p € [0, 7r/2) follows from (55). 



12 



V 

(a) 



v 

(b) 




H o 6 
(d) 



Figure 2: Moduli of kinematic (dashed), Fermi (solid), spectroscopic (dot-dashed) and astro- 
metric (dotted) relative velocities in the de Sitter universe with scale factor a(t) = e Hot . (a): 
instant comparison with respect to the function v — Hoe H ° T x', since v can be interpreted as 
a "Hubble speed" (see Remark 4.1), it is represented in dashed grey in order to compare, 
(b): retarded comparison, (c): kinematic and Fermi velocities with respect to H p. (d): 
spectroscopic and astrometric velocities with respect to HqS. 



4.2.2 Lightlike simultaneity for the de Sitter universe 



By plj ), x < X£ma X (T) = J T jf^dr = i (l-e H " T ) for any point (t,x) on the past- 
pointing horismos E~ (at proper time r of the central observer). Thus, from (53), < 



v(t, ■) < e HoT - 1 and v(-, X ) > 



By (201 



oX 



-fj— dT = y 



U(t,x) = T- — ln(l + u) 
-Wo 



So, from (41) 



S(t,x) = 



„H t 



df = 



te( T ,x) 



1 V 

H~ a l + v 



H 6 
l-H S' 



Then, by (40) 



85 



dr {l + vf 



H S(1-H S). 



Applying (39) and (58) gives, 



IKpccH 



g2HQT _ ^2H()t/> 



(l + vf-i i-(i-H sy 



e 2H r +e 2H t e ^ + ^2 + j X + ^ _ HqS ) 2 ' 



where H S e [0, 1) from (59). 



Using Remark 3.1, we can make retarded comparisons. From (54) and (58), 

1 v 

t* = t — — — In (1 + 2v) =>■ v* = , 



So, 



"kin I 



VI + V* 2 l+v ' 



"Fcrr 



l + v* 



Vy/1 + 2V 



(58) 



(59) 



(60) 



(61) 



(62) 



(63) 



Observe that all of the relative velocities for the de Sitter universe are independent of r (see 
Figure [2]). 
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4.3 Power scale factor 



In this section, we consider a power scale factor of the form a(t) = t a with < a < 1. There 
are some important particular cases, as a = 1/3, a = 1/2 (radiation-dominated universe), or 
a = 2/3 (matter-dominated universe). From ([43]) , 



(64) 



4.3.1 Spacelike simultaneity for a(t) = t a 

It follows from Q that \ < Xsmax(T) = / r i^^f=f^df = xz^C Q for any point (t,x) 



on the Fermi space slice at proper time r of the central observer, where C a 



Hence, by (64), we have that < v < v smax where, 

"smax := v (r, X smax (T)) 



1 - a 



v^r(^) 



(65) 



We have that i> sm ax is unbounded and is an increasing function of a. It reaches 1 at a = 1/3. 
% <H1>' St. W Vr^-f ^ df = and then 



l-Q 



1 1 - a 1 + a ft 



2' 2a ' 2a 'It 



2qN 



1 - a 
= C Q u, 



(66) 



where 2-Pi("> S ■> ■) is the Gauss hypergeometric function. Define the function ^(z) 
z 1 ^ 21 2 Fi (5, ; ; z 2 ) where < z < 1 and < a < 1. It is bijective and by ([66) 



withG Q («) := (i^ 1 (C Q 
From d37| and <f6Tb 



l-a v \\ 1 / a 



t s (r,x) = G a (v)r, (67) 
)) 1/Q , where the superscript _1 denotes the inverse function. 



where J a («) := C a - %^ 2 Fi (i G£ Q (u)). Setting t s (r,x) - in @ results 

in p(r, Xsmax) = C q t, the radius of the Fermi space slice at the central observer's proper time 
r, which was also found in [7J. 
By d38[) and (pi) 



f2<v 



dr = Jq,(u)t, 



(68) 



Or 



= J a (v) - (1 - a)J' a (v)v. 



(69) 



Differentiating G Q with respect to v, we obtain G„ = -\G a a ^J\ ~ G 2a , which can also be 
found by differentiating ([5]) with respect to v, taking into account that x — and '^7|. 
Integrating this differential equation and using the initial condition, G Q (0) = 1, and the 
expression of J a , gives 



J a (v) = Gi-»Vl-G^( v ) + — v 



J» = ~Gl a {v). 



(70) 



Combining (36) and (67 1 along with (70 1, (681, and (69 1, we have, 



Proposition 4.1 The kinematic and Fermi speeds of a comoving test particle relative to a 
comoving central observer in a Robertson-Walker cosmology with scale factor a(t) — t a and 
< a < 1 are given by 

t 2a 



= \ 1- 



= y/l-Gl°(v), 



(71) 
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Figure 3: Instant ((a), (b)) and retarded ((c), (d)) comparison of the moduli of kinematic 
(dashed), Fermi (solid), spectroscopic (dot-dashed) and astrometric (dotted) relative veloci- 
ties with respect to the function v — x/3r 2 / 3 ((a), (c)) and with respect to r/% 3 / 2 ((b), (d)), 
in a universe with power scale factor a(t) — t 1 / 3 . Since v can be interpreted as a "Hubble 
speed" (see Remark 4.1 1, it is represented in dashed grey in order to compare, but it coincides 
with the kinematic relative velocity in the instant comparison. 



IkFcrmiH = G 1 - a (v)y/l-Gl<>(v) + 



1 — a 



1 — a 



(1-G 2 «(v))v=t Gl a {v)v. (72) 



Remark 4.2 From (68) and since J a (v) is bijective for < v < w sma x, we have v = J a 1 (£), 
where the superscript _1 denotes the inverse function. So, by (71) and (69) (or Proposition 



4.1) we can deduce expressions for ||vkin|| and ||t>Fermi|| in terms of p/r G [0, C a ). 



Remark 4.3 An alternative expression to (72) was given in [7], 



l-a 1 



\\VFc- 



CTO 



d(J, 



(73) 



= P _ _ 

t a 2er J x a^\/a - 1 

where a = (a(r) / a(t s )) 2 = (t/£ s ) 2q = l/G 2 l Q (u). Using (66), it can be shown that v 
i rcr — ^ 1 dcr, and so (72) and (73) are equivalent. 



cr2 Q 

We turn now to examples for particular values of a. 



Example 4.1 If a = 1/3 then, from (65) we have v s 



1. From (66) we get t s ( r, \) 



[1 — u 2 ) 3 / 2 r, and by (68) we have p(r,x) = u (3 — v 2 )t. So, by (|9| and Proposition 
find 

ill =v(l + v 2 ). 



4.1 



Vki, 



V F c: 



As indicated in Remark 



4.2 



can get expressions for 



(74) 

we find v = 2 sin ^ | arcsin ( ^ J ^ , and then applying ( 74 ) we 
| and ||i>Fermi|| in terms of p/r € [0, 2) (see Figure [6]). 
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Example 4.2 For a = 1/2, i.e., for the radiation-dominated universe, (65) becomes u smax = 
tt/2. By (66), t s (T,x) — tcos 2 v, and from (68) we find p(r,x) = i v + cos i> sin i;)t. So, by 
([9]) and Proposition 4.1 we find 



smi) 



i'Fc 



i|| = (cos v + v sinii) sin v. 



(75) 



As indicated in Remark 4.2 we find that (68) defines implicitly v (p/r), and then, from ([75 
we can calculate ||vkm|| and ||«Fermi|| m terms of p/r e [0,7r/2) (see Figure[6]). 



Example 4.3 For a = 2/3, the matter-dominated universe, (65) becomes v sm ax = 2C2/3 



o V¥r(5/4) 
z r( 3 /4) 



2.62206. From (66 1 we get i s (r, x) = cd 



3 jv 



l) r, and by ( 68 ) we get p(r, x) 



■(51-1) V 1 - 4x14 (3 1 - 1) + 3) 

(|9|) and Proposition 4.1 we find 



|t'ki„|| = A/l-cd 4 (^|-l 



r, where cd (u \ m) is a Jacobi elliptic function. So, by 



II^Fc 



i !|=cd(^|-l)|| Ukin || + ^||t; kin || 2 . (76) 



As in the previous two examples, equation (68) defines implicitly the function v {p/r), and 
then, from (76) we can compute ||i>kin|| and [ppermill in terms of p/r £ [0, C2/3) (see Figure 
[61. 



4.3.2 Lightlike simultaneity for a(t) = t a 

By ([2l]), x < X^max(T) = Jo =s df = for any point (i,x) on the past-pointing horismos 
Ep (at proper time r of the central observer). Thus from (64), we have that < v < W£ max 

a 



where 
We have that 



!-7.n 



V(T, Xftna X (r)) 



1 — a 



By (20) 



is increasing with a and it reaches 1 at a = 1/2. 



(77) 



1 



dr = x 



U(T,x) = (r 1 - a -(l-a)x) i: 



1 - a 
1 v 



So, from (41 ) 



(78) 



S(t, X ) = 



tt(r,x) 



T a l + a 



1 -11- 



Now using (27), (39 1 and (40) together with (78) and (79), we have, 



(79) 



Proposition 4.2 The spectroscopic and astrometric speeds of a comoving test particle rel- 
ative to a comoving central observer in a Robertson- Walker cosmology with scale factor 
a{t) = t a and < a < 1 are given by 



Ik 



Wast 



1 — a 



-v 1 



1 -a 



(80) 



Remark 4.4 From ( 79 ) wc obtain 

1 - a 



Then, from Proposition 



4.2 



1- (l + a) 



we can express ||w S pec|| and ||vast|| in terms of S/t € 



(81) 



0. 



l+a 
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ji/2 




(b) 



(d) 



Figure 4: Instant ((a), (b)) and retarded ((c), (d)) comparison of the moduli of kinematic 
(dashed), Fermi (solid), spectroscopic (dot-dashed) and astrometric (dotted) relative veloci- 
ties with respect to the function v = X/2T 1 / 2 ((a), (c)) and with respect to t/\ 2 ((b), (d)), 
in a universe with power scale factor a(t) = t 1 ^ 2 (radiation-dominated universe). Since v 
can be interpreted as a "Hubble speed" (see Remark 4.1), it is represented in dashed grey in 
order to compare. 



Using Remark 3.1 we can make retarded comparisons. From (67) and (78 1 one can solve 



for r*. Then, with the *-version of (64) we have that v* ~ ( r .°i-a X: an d 



1 — a. 



£all = Vl " GI*(V*) ; \\v* Fe rnu\\=G 1 - a (v*)\\v* ki J + IKinl 



V. 



(82) 



Again from (64), we have r = (fx) 1 °" an< i T * = i^X) 1 " ■ So, from (67) and (78) we can 
solve v* directly in terms of v using the equation 



2F1 



1 1 — a 1 + a 



2' 2a 



2a 



w = C a , 



(83) 



where w 



Example 4.4 If a = 1/3, then from ([77 1 ve max = 1/2. From ( |78|) we get te(r,x) = (1 
2v) 3 ^ 2 t, and by (|79[), <5(t, x) = 3u (1 — v) r. So, from Proposition 



4.2 



|| V| 



spec I 



V 

1 - V 



1 1 Vast I = V (1 + «) ■ 



Taking into account Remark 



As Remark 



3.1 



from (|67j and ((7§ we have v* = ((± - l) 2 + 1 



(84) 
1/2 



4.4 



and II^Fd-mill are given by (74), with v* replacing v 



indicates, from ( 81 ) or ( 79 ) we obtain 2v = 1 — (l — 



4 J\l/2 



(84) we get ||u s: 



1+ Li!)*/ 2 and l! Vast ' 



3 t y 



and then, from 



1 



15 
3 T 



(J - l!) 1/2 , where <S/r G [0,3/4) 



3 T I 



(see Figure [6]). 

Example 4.5 For the radiation-dominated universe, with a — 



From ([78]), ^(r,%) = (1 - w) 2 r, and ([79} becomes <5(t,x) = § 



= 1/2, ([77} gives vi„ 
1 - (1 -T5*) t. So, fi 



1. 



from (80) 
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2 2.622 



(b) 



(d) 



Figure 5: Instant ((a), (b)) and retarded ((c), (d)) comparison of the moduli of kinematic 
(dashed), Fermi (solid), spectroscopic (dot-dashed) and astrometric (dotted) relative veloci- 
ties with respect to the function v = 3X/2T 1 / 3 ((a), (c)) and with respect to r /\ z ((b), (d)), 
in a universe with power scale factor a(t) — < 2 / 3 (matter-dominated universe). Since v can be 



interpreted as a "Hubble speed" (see Remark 4.1), it is represented in dashed grey in order 
to compare. 



we obtain 



Applying Remark 

< v* < v sma 



i + (i-i;) 2 



= v 1 — v 
3 



From Remark 

i-(i-H) 2/3 



3.1 

H?. So, IKJ and 



4.4 



along with ( |67[ ) and (|78j), v* is the inverse of v 

,* 

Fermi 



D* +COS V * ' 



(85) 
where 



and ( 81 ) , we get v = 1 — (1 



are given by (75), with u* replacing v. 

35^/3 

2 T> 



and then ( 85 1 results in 



i+(i-!f) 



77W3 



and 



Wast | 



= 1 II (1 



|f) 2/3 , where S/t G [0,2/3) (see Figure 



0) 



Example 4.6 For the matter-dominated universe, with a — 2/3, from (77) we have i^max = 
2. From ([78]), t e {T, X ) = (l - f) 3 r, and ([T9j) becomes 5(r,x) = § (l - (l - f) 5 ) r. So, by 



(80) we obtain 



^spec I 



pastil = 1 



V 






1 


2 


-2 


(I)'- 


5 



id)' 



(86) 



Taking into account Remark 3.1 from (67) and (78) we obtain v* as the inverse of v = 
„.+ 2c d(^|-i) ' where 0<v * < ^max = 2C 2/3 ps 2.62206. So, ||u*J| and K ermi || are given 
by (76), consideri ng v * instead of v. 
From Remark 



5S\ 4 /5 



4.4 



and (81) we get ^ = 



i-(i-H) 
i+(i-fl) 4/5 

Figure M. 



and 



= l-|i-2(l 



1 - - II) and then © § iveS ll«-po 



5 A^ 3 / 5 
3 



_ 5S\ 
3 t) 



1 - |f) 4/5 , where S/t e [0,3/5) (see 
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Figure 6: Moduli of the relative velocities in a universe with power scale factor a(t) = t a with 
a = 1/3 (left), a — 1/2 (middle) and a — 2/3 (right). Top: kinematic (dashed) and Fermi 
(solid) velocities with respect to p/r. Bottom: spectroscopic (dot-dashed) and astrometric 
(dotted) velocities with respect to S/t. 



5 Functional relationships for relative velocities 

In this section, we allow the radial motion of test particles to be non geodesic. Throughout, we 
assume that the scale factor a(t) is smooth and increasing. We employ both Fermi coordinates 
and optical coordinates to derive functional relationships between the Fermi, kinematic, 
astrometric, and spectroscopic relative velocities of test particles. First, we relate Fermi 
and kinematic relative velocities to each other and spectroscopic and astrometric relative 
velocities to each other, and consider some geometric consequences of these relationships. 
Then, in Subsections |5 . 1 1 and |5 . 2| we use Fermi coordinates to find retarded relationships for 
the spectroscopic and astrometric relative velocities respectively. 

We begin with a relationship between the Fermi and kinematic relative velocities that 



follows immediately from ( 14 1 



VFc 



= V-9tt(t, p) Wkin- (87) 



Analogously, in the context of lightlike simultaneity, combining (30) with (27) gives 

1+9tt(t, 6) ± 1 1 V ast | 



1 - 9tt(t, S) =F 1 1 «ast 



\S P , (88) 



where in the case that dS/dr > (i.e., in the case of a receding test particle), the positive 
sign in the numerator and negative sign in the denominator are chosen, and the opposite 
choices of signs are taken when d5/dr < (i.e., in the case of an approaching test particle). 



From (87) and (88) we may formulate the following proposition. 



Proposition 5.1 For a Robertson-Walker spacetime with scale factor a(t) that is a smooth, 
increasing function oft, the kinematic and Fermi speeds of any test particle undergoing radial 
motion with respect to a comoving observer determine the Fermi metric tensor element g TT 
at the spacetime point of the particle, via 

/ \ II ^Fermi 1 1 

9tt{t,P) = n — , 

||«kin|| 

provided the denominator is not zero. Similarly, a measurement of the astrometric and 
spectroscopic speeds of a receding test particle relative to the comoving observer determine 
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the metric tensor element g TT in optical coordinates at the spacetime point of the particle via 

- i X s 1 - bast || || || 1 + ll^ast || 

g TT (T,d) = —— j| Kpcc - T— jj M- 

l+||Vspec|| 1 + Fspecll 

Remark 5.1 It follows from Proposition |5.1| that in principle, knowledge of either pair of 
the relative velocities for moving test particles at each spacetime point uniquely determines 
the geometry of the two dimensional spacetime via Q and p6| ) , and therefore the scale factor 
a(t). Since the affine distance (i.e the optical coordinate 5) can be measured by parallax, 
and the frequency ratio can be found by spectroscopic measurements, the astrometric and 
spectroscopic relative velocities can be determined solely by physical measurements, and so, 
they could confirm or contradict assumptions about the value of a(t) for the actual universe. 



Proposition 5.1 along with results of the preceding section allow us to calculate g TT (T,p) 



in terms of the curvature-normalized coordinates. For example, equations (67), (71) and 
Proposition |4. 1| immediately give the following result: 



Corollary 5.1 For a Robertson- Walker spacetime with scale factor a(t) = t a , with < a < 
1, the metric tensor element g TT is given by 

9rr (r, p) = - (a 1 -* (v)+ 1 ^ V^-GlHv) w) 2 , 

where v — v(t,x) ~ a xl rl ~ a , an d X * s ihe unique coordinate such that (G a (v)T,x) in 
curvature-normalized coordinates represents the same event as (r, p) in Fermi coordinates. 



As an illustration of Corollary |5.1| for the radiation-dominated universe, i.e., when a = 
1/2, we immediately have from Example 4.2 that g T r{T,p) = — (cosw + vsinv) 2 , and it may 



be verified using the coordinate transformation formulas in [7] that this expression is the 



same as (10). 



5.1 Retarded comparison for the spectroscopic relative velocity 

The redshift of a radially moving test particle relative to the central observer is determined 



by the frequency ratio (28). In Fermi coordinates, the 4-momentum tangent vector field for 



the photon is defined on the light ray A and it is given by 



P = i 



d_ 



. d 

P dp = 



V d v d 

Ot dp ' 



(89) 



where the overdot represents differentiation with respect to the affine parameter of A, and, 
analogous to £ in ( 11 ), V := sJ—g TT f is the energy of the photon as measured by a stationary 



observer (in Fermi coordinates), i.e. an observer with 4- velocity 



Combining fl28h, (89) and (11), gives 



V-9tt <9t- 



Hit) 



S* ± v 7 ^* 2 - i) = — (s* ± \/£* 2 - i 



(90) 



where the + or — sign is taken for a receding or approaching (in Fermi coordinates) test 
particle respectively, the energy £* is calculated for the test particle located at qi = (r*,p*), 
and Vq/p is the frequency ratio of a photon emitted from a stationary observer (whose 4- 
velocity is given by (12)) from the event qg. Using (15), (90) may be expressed as 



^£*(l±\\v } 



kin I 



(91) 
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where v^- ln is the kinematic relative velocity of u' e measured from the event p* of j3 with Fermi 
coordinates (r* , 0) (see ( 14 ) ) . The term Vq/v may be calculated by considering a test particle 
comoving with the Hubble flow with fixed curvature-normalized coordinate Xi f° r which the 
left side of (91) is known. In that case, the left side of (91) is given by (42), i.e., 



v 
v 



o(r)_ 
a(tt)'' 



(92) 



where tt = ttij, x) is determined implicitly by (20), and the coordinates (tg,x) an d (t* , p*) 
describe the same event qg. The 4- velocity of the comoving test particle is given by 



01, 



Op 



Qi 



. Thus, 



y/-grr(T*, P *) 9l 



d 




11 





V-9tt(t*,p*)^t(U, X ). (93) 



The partial derivative on the right side of (93) may be calculated as follows. Since \ is fixed, 
we may write ti — te(r). Starting with the coordinate transformation formula, t(£^(t),x) = 
t* , we have by the chain rule, 



dr* 



dr dt 
dT^Tr 



dr 

dtp. 



dr* dr 
dr dtp 



(94) 



Next we find dr* /dr. Since the vector field P is tangent to A, using (89) we find 



df = 



1 



.dp 



V-(Jtt{t(p),p) 
which determines r as a function of r* with 



p*(r*) 



1 



dr 
dr* 



= 1 



1 



dp* 



y/-grr{T*, P *)dT* 



o V-9tt(t(p),p) 



= i + Ktall >o. 



dp, (95) 



(96) 



It follows that (95) determines the inverse function, r*(r) as well. On the other hand, 

ain 

(97) 



differentiating (20) with respect to r and taking into account that \ is constant, we obtain 

a{tt) 



dt e 
~dr 



o(t) 



Using (96) and (97) to calculate the right side of (94), we have jsr(tt,x) — \7t5i * \\\ ■ 
Thus, from (93), we have for the comoving test particle, 



£* 



a ( T )V-9Tr{T*,p*) 

<u) (i + KJI) ' 



(98) 



Combining ( 98 1 and ( 92 ) with ( |9l[ ) , yields the frequency ratio of a photon emitted by a 
stationary observer (in Fermi coordinates) at qi — (r* , p* ) and received by the Fermi observer 
at p = (t,0), 

^ = , 1 (99) 



Now combining (99) with (91) gives the following expression for a radially moving particle 

(100) 



(not necessarily comoving with the Hubble flow): 

£*(l±Knl 



v V~9tt(t*,p*)' 
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The corresponding spectroscopic relative velocity may be computed and expressed in terms 
of the kinematic relative velocity for a radially moving particle, directly from Equation 
(12)] and (91 1 above, as 

<"' / "» 2 - 1 5„ (101) 



^spcc 



(i//i/) 2 + 1 p 



Specializing, for the sake of simplicity, to the case of a radially receding test particle ( 101 ) 
may be rewritten as 



v 
v 




Now, combining (102 1 with (100 1 yields an expression for ||t£ in || in terms of ||v S p e 



\/-9tt(t*,p*) l+\\v spec \\ 



(102) 



(103) 



5.2 Retarded comparison for the astrometric relative velocity 

The modulus of the astrometric relative velocity of a radially moving test particle relative 
to /3 is \dS/dr\, where S is the affine distance from p to qe (see Figure [I]). To compute this 
we use the fact that A(<5) = qt, where A(£) is an affinely parameterized null geodesic from 



A(0) =pto 



with ^ 



— I 
dr \p 



0,, 



(see [3 Propositions 6 and 7]) and whose 



tangent vector field P is given by (89). Therefore, p* = J Q pd£, where the overdot represents 
differentiation with respect to £, and thus 



dp*_ 



P(S) = -V-9rr(r*, P *) t(6) = g(u ,P qe ) = -V(q e ), 



where Uq is the stationary observer at qi . Note that the previous expression is positive since 
the photon is backward traveling. In order to calculate V(qe), consider a forward moving 
photon traveling from qi = (r*,p*) to p — (r, 0) following the future-pointing null geodesic 

'o - i 



A(f) = X(8 - 0- Then (28) and (99) give -V{q t ) 



and thus, 



dS_ 
~dp* 



= V-9tt{t*,p*). 



(104) 



Using the chain rule together with (96), (104) and the *-version of ([9]), it follows that the 
modulus of the astrometric relative velocity of u' g observed by u is given by 



Fast = 



dS 



dS 
Jp* 



dp* 



dr* 



dr* 
dr 



dp* 



dr* 



\J-9tt{t*,P*) 



^-g T Ar',P') dT * 

which may be expressed in terms of the Fermi or kinematic relative velocity of u' £ measured 
from the event p* of j3 with Fermi coordinates (t* , 0) as, 



I Vast || = 



^-9rr(r*,p*)\\v* FcImi \\ _ -<? T r(T*,P*)K in || 



1± 



1 ± \\v 



(105) 
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where the positive sign is used for a radially receding particle, and the negative sign for a 



radially approaching particle, and where for the second equality we used (87) 



Finally, we can find relationships between the astrometric and spectroscopic relative veloc- 
ities. For ease of exposition, we specialize to the case of a receding test particle. Rearranging 



( 105 1, we have 



Kinll = 



I Wast | 



-(r*,p*) - 1 1 ^ 1 1 ' 
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and then from (1031, wc find 



Kurt 1 1 = ~9tt{t*,p*) 1 



£* 



\ V-9rr(T*,p*) V l + |Kpec|| J ' 

This may be compared with the following expression in optical coordinates, which follows 
from ([27| and Q (see also (88])), 



I Wast | 



= ~o 3rr(T,(5) + 



l- Ik 



spec I 



1 + IkspecH 



6 Concluding remarks 

We have found general expressions for the Fermi, kinematic, astrometric, and spectroscopic 
velocities of test particles experiencing radial motion relative to an observer comoving with 
the Hubble flow (called the central observer) in any expanding Robertson- Walker cosmology. 
Specific numerical calculations and formulas were given for the de Sitter universe and cos- 
mologies for which the scale factor, a(t) = t a with < a < 1, including, in particular, the 
radiation-dominated and matter-dominated universes. 

In Proposition |5.1| we showed how pairs of these relative velocities, for arbitrary test 
particles, determine the leading metric tensor coefficients in Fermi or optical coordinates, 
and in Remark |5.1| we observed that the astrometric and spectroscopic relative velocities 
could be used, at least in principle, to measure the scale factor a(t) of the actual universe. 
An analog of Hubble's law for the astrometric relative velocity of comoving particles was 
given in Proposition |3.1| 

Of the four relative velocities, only the Fermi relative velocity of a radially receding test 
particle can exceed the speed of light, and this is possible at a spacetime point (t, p), in Fermi 
coordinates, if and only if —g T T{T,p) > 1. Examples of superluminal Fermi velocities were 
given in Section [4j 

Under general conditions, the Hubble velocity of comoving test particles also become 
superluminal at large values of the radial parameter, x, and this is taken as a criterion 
for the expansion of space in cosmological models, and for the actual universe. By way 
of comparison, the Fermi relative velocity has both advantages and disadvantages to the 
Hubble velocity. For comoving particles, both velocities measure the rate of change of proper 
distance away from the observer with respect to the proper time of the observer. But for 
the Fermi velocity, the proper distance is measured along spacelike geodesies, while for the 
Hubble velocity the proper distance is measured along non geodesic paths. In this respect 
the Fermi velocity is more natural and more closely tied to the observer's natural frame of 
reference, i.e., to Fermi coordinates in which locally the metric is Minkowskian to first order 
in the coordinates. In addition, the notion of Fermi relative velocity (as well as the other 
three relative velocities considered herein) is geometric and applies to any spacetime, while 
the Hubble velocity is specific to Robertson- Walker cosmologies. 

In its favor, the Hubble velocity is defined at all spacetime points, whereas the Fermi 
velocity makes sense only on the Fermi chart of the central observer. For the case of non 
inflationary cosmologies, such as when the scale factor follows a power law, a(t) = t a , the 
Fermi chart is global [7], and the Fermi velocity does not suffer this disadvantage. 

However, for inflationary cosmologies such as the de Sitter universe, the Fermi chart is 
valid only up to the cosmological horizon 11511161 17]. For this particular example, the behav- 
ior of the Fermi and Hubble velocities of comoving particles is qualitatively different. Fermi 
relative velocities do not reach half the speed of light, while Hubble velocities become super- 
luminal, as shown in Figure [2] This is because the notions of simultaneity are different for 
the corresponding coordinate systems, and the two velocities are "measuring" the expansion 
of different hyperspaces. It would be interesting and possibly useful to have a rigorous and 
universal mathematical criterion for expansion of space. 
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